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Abstract 

We have simultaneously measured the energy spectra of neutrons and protons 
emitted in the non-mesonic weak decays of ^He and ^ 2 C hypernuclei produced via 
the {n + ,K + ) reaction with much higher statistics over those of previous experi- 
ments. The neutron-to-proton yield ratios for both hypernuclei at a high energy 
threshold (60 MeV) were approximately equal to two, which suggests that the ratio 
of the neutron- and proton-induced decay channels, T n (An — ► nn)/T p (Ap — > np), is 
about 0.5. In the neutron energy spectra, we found that the yield of the low-energy 
component is unexpectedly large, even for ^He. 
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1 Introduction 



In free space, a A particle decays dominantly associated with a pion in the 
final state as A —>■ Nir. In the case of A bound in a nucleus, a A hypernucleus, 
A is not only possible to decay as in the free space (mesonic weak decay), but 
is also able to stimulate a nucleon in the nucleus as A "JV" — » nN (non- mesonic 
weak decay, NMWD). The NMWD gives a unique opportunity to study the 
weak interaction between baryons, because this strangeness non-conserving 
process is purely attributed to the weak interaction. 

In the NMWD of a A hypernucleus, there are two decay channels, Ap — > np 
(T p ) and An — >• nn (T n ). The ratio of those decay widths, T n /T p , is an impor- 
tant observable used to study the isospin structure of the NMWD mechanism. 
For the past 40 years, there has been a longstanding puzzle that the exper- 
imental T n /T p ratio disagrees with that of theoretical calculations based on 
the most natural and simplest model, the One-Pion Exchange model (OPE). 
In this model, the AN — > nN reaction is expressed as a pion absorption pro- 
cess after the A — > Nn decay inside the nucleus. Since the OPE process is 
tensor-dominant and the tenser transition of the initial AN pair in the s-state 
requires the final nN pair to have isospin zero, the T n /T p ratio in OPE process 
becomes close to 0. However, previous experimental results have indicated a 
large T n /T p ratio (~1) [1,2]. 

This large discrepancy between the OPE-model predictions and the experi- 
mental results has stimulated many theoretical studies: the heavy meson ex- 
change model [3,4,5], the Direct Quark model [6,7,8,9,10] and the two-nucleon 
(27V) induced model (ANN -> nNN) [11]. After K. Sasaki et al. pointed out 
an error about the sign of the Kaon exchange amplitudes in 2000 [10], those 
theoretical values of the r n /r p ratio have increased to the level of 0.4~0.7 
[12]. 
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On the other hand, the experimental data still have large errors {T n /T p = 
0.93 ± 0.55 for ^He [1]), and it is hard to draw a definite conclusion on the 
T n /T p ratio. So far, several nucleon energy spectra from hypernuclear decay 
have been reported. Most of the experiments measured only the proton energy 
spectra, because of the difficulty in detecting neutrons. The T n /T p ratio has 
been estimated through a comparison of the measured proton spectrum with 
that of the intra- nuclear cascade calculation by changing the T n /T p ratio [13]. 
In this method, there are large uncertainties due to the indirect evaluation of 
neutrons. 

Recently, the neutron energy spectra from ^ 2 C and ^ 9 Y were measured with 
high statistics [14]. For ^ 2 C, the T n /T p ratio was extracted from the neuron- 
to-proton yield ratio, N n /N p , estimated by the neutron energy spectrum com- 
bined with a proton energy spectrum measured in another experiment [13]. 
The energy loss of protons in the target, however, was not able to be cor- 
rected for. Since there were ambiguities coming from different detection energy 
thresholds in the two measurements, the value might be affected by a large sys- 
tematic error. We therefore simultaneously measured the neutron and proton 
spectra from the decay of A hypernuclei with a proton energy loss correction, 
and with much higher statistics over those of previous experiments. 

When we compare the "measured" T n /T p ratio with that obtained in theoret- 
ical calculations, the most serious technical problem was a treatment of the 
re-scattering effect in the residual nucleus, the so-called Final State Interac- 
tion (FSI). The nucleon energy spectrum emitted from the two-body NMWD 
process, AN — > nN, must originally have a broad peak at half of the Q-value 
(~ 76 MeV). The FSI effect distorts the nucleon spectra enhancing the low 
energy region, in which the strength of the effect depends on the mass number 
of the hypernuclei. In order to minimize the FSI effect, we selected a light s- 
shell hypernucleus, ^He. In the s-shell hypernucleus, initial relative AN states 
must be S states, whereas in a p-shell hypernucleus it is possible to be P 
states. According to Ref. [15], the transitions coming from the P states were 
predicted to be small, though we need experimental confirmation for that. In 
order to investigate the p-wave effect, we also performed the same experiment 
for a typical light p-shell hypernucleus, ]fC In this paper, we show those 
energy spectra for both hypernuclei and present the neuron-to-proton yield 
ratios, N n /N p , with a high energy threshold (60 MeV) in order to get rid of 
the contribution from the FSI effect. 

On the other hand, the nucleon yields in the low energy region should be 
sensitive to the strength of the FSI effect. In the proton measurement there 
is a cutoff energy, typically at 30 ~ 40 MeV, due to the energy loss inside 
the target and the detector, whereas in the neutron measurement we have 
experienced sensitivity in the whole energy region. In this paper, we also show 
the mass number dependence of the "neutron" energy spectra for A = 5 and 12, 



with that for A = 89 [14], which provides valuable information to investigate 
the FSI effect. The possible existence of a multi-nucleon induced process in 
which the Q-value can be distributed over three or more nucleons in the final 
state has been discussed theoretically (such as 2iV-induced process), though 
there has been no experimental evidence so far; the nucleons emitted from 
this process should also distribute in the low energy region. Thus, the whole 
shape of the neutron energy spectrum also brings useful knowledge on the 
multi-nucleon induced NMWD mechanism. 



2 Experimental method 



The present experiments (KEK-PS E462/E508) were carried out at the 12- 
GeV proton synchrotron (PS) in the High Energy Accelerator Research Orga- 
nization (KEK). Hypernuclei, ^He and ^ 2 C, were produced via the (ti + ,K + ) 
reaction at 1.05 GeV/c on 6 Li and 12 C targets, respectively. Since the ground 
state (g.s.) of ^Li is above the threshold of ^He + p, the g.s. promptly de- 
cays into 5 A He emitting a low-energy proton. The 6 Li (tt + ,K + ) ^Li reaction 
was therefore employed to produce ^He. The hypernuclear mass spectra were 
calculated by reconstructing the momenta of incoming ir + and outgoing K + 
using a beam-line spectrometer composed of the QQDQQ system and the 
superconducting kaon spectrometer (SKS) [16], respectively. 



Particles emitted from the decays of A hypernuclei were detected by the decay- 
particle detection system, as shown in Fig. 1. It was composed of plastic scin- 
tillation counters (Tl, T2, T3 and T4) and a multi-wire drift chamber (PDC). 
The neutral particles, neutrons and gammas, were identified by means of the 
time-of-flight (TOF) technique between Tl and T4; in addition, T3 was used 
as a veto counter in order to reject charged particles. The charged particles, 
protons and charged pions, were identified by utilizing dE/dx on T2, total 
energy deposit on sequentially fired counters (T2, T3 and T4), and the TOF 
between T2 and T3. Their trajectories were measured by PDC. 




£+ to S^ s 



Drift chamber (PDC) 




40cW 



Fig. 1. Schematic view of the decay-particle detection system. 



3 Analysis 



Figure 2 shows the inclusive excitation-energy spectra of ^Li and ^ C. The g.s. 
peak is clearly seen in our excitation spectra for both hypernuclei. We selected 
events from the hypernuclear production by gating the g.s. peak, as shown in 
the figure. 



The yields of ^He and A C are, respectively, about 4.6 x 10 and 6.2 x 10 , 
which were one order-of-magnitude higher than those of previous experiments. 
There is a constant background below the g.s. peak; the background for ^He 
and ]^C within the g.s. gate were estimated to be as small as 8% and 3%, 
respectively. We estimated the contribution from this constant background in 
the nucleon spectra and subtracted it from the g.s. gated nucleon spectra. 
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Fig. 2. Inclusive excitation-energy spectra of (a) ^Li and (b) ^ 2 C. 

In the neutral decay particle analysis, the TOF between the reaction vertex 
and the T4 counter, TOF n , is expressed as TOF n = Tta — Tri — TOF\, — At, 
where Tta and Tt\ are the timings of T4 and Tl (time zero), respectively; 
TOFb means the flight time of the beam pion; At corresponds to the time 
delay due to the hypernuclear lifetime. Since At cannot be determined for each 
event individually, we replaced At with the mean lifetime of the hypernuclei 
to avoid a systematic energy shift in the neutron energy spectra. Here, we used 
the lifetimes of r = 278±]q P s for 5i He and r = 212^ ps for 1;i 
through the lifetime analysis of the present experiment [17]. 
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The inverse of the neutron velocity, 1/(3, was calculated as 1/(3 = TOF n / 
L x c, where L is the flight length and c is the light velocity. The top of 
Fig. 3 shows the measured 1//3 spectra with a light-output threshold of 2 
MeVee (MeV electron equivalent) for ^ 2 C. The energy scale is displayed on 
the upper horizontal axis. The gate region of the neutron was defined as 1.97 
< 1/(3 < 9.73, which corresponds to 5 < E n < 150 MeV of neutron energy. The 
spectrum shows good 7/neutron separation. The extremely low background 
level of the 1/(3 < region in the spectrum indicates that the accidental 
background is negligibly small. The neutron energy resolutions for ^He and 
^ 2 C were estimated from the width of the 7-ray peak. They were, respectively, 



about 10 MeV and 7 MeV in full width at half maximum (FWHM) at E„ 
MeV. 



=75 



The bottom of Fig. 3 shows the particle identification (PID) of the charged 
particles emitted from the decay of ^He. This plot was made by the PID 
function derived from the dE/dx, the total energy and the particle velocity. 



The charged particles (charged pions, protons and deuterons) were separated 
well. It is noteworthy that deuterons were separated from the protons for the 
first time in such counter experiments for the NMWD. The proton peak was 
fitted by a Gaussian function, and the proton gate was set at the ± 2 a region 
of the Gaussian, as shown in the figure. The loss of protons (4.6%) by this tight 
cut was taken into account in the analysis efficiency The proton kinetic energy 
was determined by the TOF between T2 and T3, which must be affected by 
the energy losses inside the target and in the T2 counter. Those energy losses 
were estimated event-by-event from the path length inside the target and the 
T2 counter, which was calculated from the reaction vertex and the trajectory. 
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Fig. 3. Particle identifications for neutral (top) and charged (bottom) decay par- 
ticles. The top shows the 1/(3 spectra of neutral particles emitted from the decay 
of ^ 2 C with 2 MeVee threshold. The bottom shows the PID function of charged 
particles emitted from the decay of 



*He. 



The numbers of the neutron and proton per NMWD at energy E, N n (E) and 
N p (E), are expressed as 



N n , p (E) 
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where Yuy is the number of g.s. formations in the inclusive spectrum; Y n (E) 
and Y p (E) are the numbers of detected neutrons and protons from the decay. 



The acceptances for the neutron and proton detection system, Q n and fl p , 
were estimated respectively to be 26.6 ± 0.8% and 10.0 ± 0.3% by a Monte- 



Carlo simulation. Only for the proton analysis, we selected the central parts 
of the segmented T2 and T3 counters of the top and bottom coincidence arms 
in order to suppress ambiguity of the acceptance estimation due to escape 
from the side edge of the T4 counter. The difference between Q n and Q p is 
attributed to this selection. 

The non-mesonic decay branching ratios of ^He and ^ 2 C, b nm , were evaluated 
as 0.429 ± 0.012(stat) ± 0.005(sys) and 0.768 ± 0.005(stat) ± 0.012(sys), 
respectively, by n~ and it branching-ratio analyses of the present experiment 

[17,18]. 

The energy-dependent detection efficiency for the neutron, e n (E), was esti- 
mated by the Monte-Carlo simulation code DEMONS [19], which is based on 
CECIL [20] and is applicable to multi-element neutron detectors. We com- 
pared the calculated efficiencies with those of various existing experimental 
data [14], and found that the systematic error does not exceed 6% for in- 
tegrated yields above 10, 20, 30 and 40 MeV. According to a Monte-Carlo 
simulation for the proton, there is no experimental sensitivity to a proton be- 
low about 30 MeV, due to the energy loss, and the detection efficiency was 
gradually reduced from 50 MeV to 30 MeV. The energy-dependent detection 
efficiency for the proton, e p (E), was calculated using the simulation. 

There is a non-negligible nucleon background due to the pion absorption pro- 
cess in which 7r _ 's from the mesonic decay of A hypernucleus are absorbed by 
the materials around the target. The background was estimated by assuming 
that the shape of the nucleon spectra from this 7r" absorption process is the 
same as that from the ir~ decay of A (A — > 7i~p) formed via the quasi-free 
formation process. Here, we regarded the 7r~ branching ratio of the quasi-free 
A decay as that of A in the free space, and estimated the contribution using 
the 7T" branching ratios, 0.359 ± 0.009 for ^He [17] and 0.099 ± 0.011 (stat) 
± 0.004(sys) for fC [13]. 



4 Results and Discussion 



The top and middle of Fig. 4 show the neutron/proton energy spectra from 
the decay of ^He and ^ 2 C, respectively. They were normalized per NMWD. 
It seems that the neutron spectra for both hypernuclei have a similar shape 
to those of protons above the proton energy threshold of 30 MeV, and the 
neutron yields were about twice higher than those of protons. Thanks to the 
event-by-event energy-loss correction for proton, we can directly compare the 
neutron energy spectra with those of the protons at the same energy threshold. 
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Fig. 4. Neutron (filled circle) and proton (open circle) energy spectra per NMWD of 
(a) ^He and (b) ]\C, with the neutron spectrum of (c) ^ 9 Y obtained in the previous 
experiment [14]. The errors are statistical. The dashed histogram in the top figure 
shows the neutron spectrum per NMWD of ^He, calculated by Garbarino et al, in 
which the FSI effect and the 2./V-induced process were taken into account [22]. 



It should be noted that the partial decay rate of the ^He — > n+a process, which 
emits a monochromatic high-energy neutron (~ 135 MeV), is not negligible. 
The decay rate was reported as 0.049 ± 0.01 T w - [21], where r„- is the n~ 
decay width of ^He. Though it is hard to observe the corresponding peak due 
to the limited resolution, ~ 25 MeV (FWHM) at such a high-energy region, 
one can see an appreciable neutron yield at the high-energy region (120 ~ 150 
MeV) in the neutron spectra of ^He, which can be considered to be the events 
from this process. 

The FSI process and the possible multi-nucleon induced process (such as the 
2iV-induced process) tend to enhance the low-energy region in the nucleon 
energy spectra. In order to reduce the effects from those processes, we set the 
energy threshold to be as high as 60 MeV, which is near to half of the Q-value 
(~ 76 MeV) of the AN — > nN decay process. The ratios of the yields between 
neutrons and protons, N n /N p , for ^He and ^ 2 C above 60 MeV, were obtained 
respectively as follows: 



N n /N p QHe)=2.17±0.15(stat) ±0.16(sys) (60 < E < HOMeV), (2) 
N n /N p (%C) = 2.00 ± 0.09(stat) ± 0.14(sys) (E > 60MeV), (3) 

where we set the upper energy limit of 110 MeV only for ^He in order to 
avoid the contamination from the ^He —>■ n + a process. Systematic errors are 
predominantly coming from an ambiguity of the absolute neutron detection 
efficiency (~6%). Due to the similar spectrum shape for neutrons and protons, 
the N n /Np ratios of both hypernuclei are quite stable for a change of the 
energy thresholds. Recently, an independent experimental result of the N n /N p 
on ^ 2 C was reported as 1.73±0.22 [14], which agrees with our result within the 
errors. If the FSI effect and the contribution of the 2iV-induced process can 
be neglected, the relation between the T n /T p and the N n /N p for AN — > nN 
decay process is approximately expressed as 

N n /N p = 2 x (ryr p ) + 1. (4) 

When we applied the above relation to the obtained N n /N p , the T n /T p ratios 
for both hypernuclei were obtained as about 0.5~0.6. The result excludes 
earlier experimental results that the T n /T p is close to unity [1,2]. It also rules 
out a theoretical calculation based on the OPE model, which predicts that the 
value should be as small as 0.1. Our result supports recent calculations based 
on short-range interactions, such as the heavy-meson exchange and the direct 
quark exchange. Moreover, we could not observe significant difference of the 
T n /T p ratio between the s-shell and p-shell hypernuclei. It is suggested small 
contribution from the p-wave in this process. 

The bottom of Fig. 4 shows the neutron spectrum per NMWD of ^ 9 Y ob- 
tained in a previous experiment [14]. With this spectrum, the mass- number 
dependence of the neutron energy spectra for A = 5, 12 and 89 is shown. The 
integrated numbers of neutrons per NMWD of ^He, ^ 2 C and ^ 9 Y are listed 
in Table 1 for energy thresholds of 10, 20, 30 and 40 MeV. The high-energy 
region is suppressed with the increase of mass number while enhancing the 
low-energy component. This tendency can be naturally interpreted as being 
the FSI effect. 

If the liV-induced process (AN — ► nN) dominates the NMWD and the FSI 
effect is negligible, the neutron energy spectrum should have a broadened peak, 
due to the Fermi motion, at about one half of the Q- value of the process (~75 
MeV). We simply expected that this would be the case for ^He due to its small 
mass number. Garbarino et al. reported a calculated nucleon energy spectrum 
for ^He with the FSI effect and assuming that r n /r p = 0.46 and T2n/^in — 
0.20, where r 1A r and r 2 7v denote the decay widths of the IN- and 2iV-induced 
processes [22]. It is overlaid on the top of Fig. 4 (dashed histogram) on the same 
scale. The observed spectrum was significantly different from the calculated 
energy spectrum. The calculated spectrum still has its maximum at about 
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Table 1 

Total number of neutrons per NMWD of ^He, A 2 C and ^ 9 Y [14] with energy thresh- 
olds of 10, 20, 30 and 40 MeV. 





Total number of neutrons per non-mesonic decay 




*He 


A ^ 


»Y [14] 


E n 


> 10 MeV 


1.398±0.052±0.096 


1.266±0.020±0.090 


1.36±0.08±0.09 


E n 


> 20 MeV 


1.168±0.044±0.080 


1.024±0.018±0.072 


0.89±0.06±0.06 


E n 


> 30 MeV 


1.036±0.039±0.071 


0.823±0.016±0.058 


0.67±0.06±0.04 


E n 


> 40 MeV 


0.918±0.036±0.063 


0.687±0.015±0.048 


0.58±0.05±0.04 



Q/2, whereas the experimental spectrum levels off in the region from 20 to 80 

MeV. 

It should be noticed that there are so far two types of the 27V-induced NMWD 
models. One is a three-body reaction in which the Q- value can be distributed 
over three nucleons in the final state. The other is based on the correlated two- 
nucleon absorption of a virtual pion emitted from the weak vertex of ANtt, 
in which the energy of the emitted nucleon is small, since the virtual pion is 
close to on-shell [22] . Since most of the Q- value in this 2iV-induced process is 
distributed to two nucleons, the shape of the energy spectra should be close to 
that in the liV-induced process (AN — > nN). Thus, the latter model adopted 
in the calculated spectra has a much smaller effect to enhance the low-energy 
region than the former model. This suggests the importance of the former 
2iV-induced model. Another speculation could be that the FSI effect is much 
stronger than that in the calculation by Garbarino et al. 



5 Conclusion 



We measured the energy spectra of neutrons and protons emitted from the 
non-mesonic decays of ^He and \ 2 C with high statistics. Both the neutron and 
proton spectra show a similar shape. The obtained neutron-to-proton yield 
ratios, N n /N p , for both hypernuclei are approximately equal to two with high 
energy thresholds, which are consistent with a simple estimation assuming 
r„/r p ~ 0.5. This result rules out a theoretical calculation based on the OPE 
model, and supports recent calculations based on short-range interactions. 
In the present experiment, the high statistics of the observed nucleons has 
enabled us to carry out a coincidence analysis of the two nucleons, n + n and 
n + p, which is now in progress. The analysis will provide a definitive T n /T p 
ratio. In the obtained neutron spectrum for ^He, the spectral shape was not 
consistent with the simple expectation that there is a broad peak at one half 
of the Q-value. This indicates the importance of the multi-nucleon induced 
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process in the NMWD or/and a large FSI effect, even for ^He. This result 
calls for more detailed studies concerning the 2iV-induced NMWD process 
and the FSI effect. 
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